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, having a wide range of disk-to-central object mass ratios. Global 
eigenmodes with four distinctly different characters have been 
; identified using numerical, nonlinear hydrodynamic techniques. 
The mode that appears most likely to arise in normal star formation 
settings, however, resembles the “eccentric instability” that has 
been identified earlier in thin, nearly Keplerian disks: It presents an 
open, one-armed spiral pattern that sweeps continuously in a trail- 
- ing direction through more than 2-tt radians, smoothly connecting 
the inner and outer edges of the disk, and requires cooperative 
motion of the point mass for effective amplification. This particular 
: instability promotes the development of a single, self-gravitating 
■ clump of material in orbit about the point mass, so its routine 
appearance in our simulations supports the conjecture that the 
_ eccentric instability provides a primary route to the formation of 
short-period binaries in protostellar systems. 
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THREE-DIMENSIONAL RADIATIVE TRANSFER CAL- 
CULATIONS ON AN SIMD MACHINE APPLIED TO 
ACCRETION DISKS. H. Vath, Department of Physics and 
Astronomy, Louisiana State University, Baton Rouge LA 70803, 
USA. 

1 We have developed a tool to solve the radiative transfer equation 
ifor a three-dimensional astrophysical object on the SIMD computer 
"MasPar MP- 1. With this tool we can rapidly calculate the image of 
Isuch an object as seen from an arbitrary direction and at an arbitrary 
1 wavelength. Such images and spectra can then be used to directly 
'compare observations with the model. This tool can be applied to 
- many different areas in astrophysics, e.g., HI disks of galaxies and 
_ polarized radiative transfer of accretion columns onto white dwarfs. 
Here we use this tool to calculate the image and spectrum of a simple 
model of an accretion disk. 
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DYNAMICS OF FLUX TUBES IN ACCRETION DISKS. 
E.T. Vishniac and R. C. Duncan, Department of Astronomy, The 
University of Texas, Austin TX 78712, USA. 

| The study of magnetized plasmas in astrophysics is complicated 
f by a number of factors, not the least of which is that in considering 
! magnetic fields in stars or accretion disks, we are considering 
-plasmas with densities well above those we can study in the labo- 
ratory. In particular, whereas laboratory plasmas are dominated by 
the confining magnetic field pressure, stars, and probably accretion 
disks, have magnetic fields whose p (ratio of gas pressure to mag- 
netic field pressure) is much greater than 1 . Observations of the Sun 
Suggest that under such circumstances the magnetic field breaks 
apart into discrete flux tubes with a small filling factor. On the other 
hand, theoretical treatments of MHD turbulence in high-P plasmas 
jend to assume that the field is more or less homogeneously distrib- 
uted throughout the plasma [11. 

Here we consider a simple model for the distribution of magnetic 
'flux tubes in a turbulent medium. We discuss the mechanism by 


which small inhomogeneities evolve into discrete flux tubes and the 
; size and distribution of such flux tubes. We then apply the model to 
' accretion disks. We find that the fibrilation of the magnetic field 
f does not enhance magnetic buoyancy. We also note that the evolu- 
: tion of an initially diffuse field in a turbulent medium, e.g., any 
: uniform field in a shearing flow, will initially show exponential 
growth as the flux tubes form. This growth saturates when the flux 
tube formation is complete and cannot be used as the basis for a self- 
1 sustaining dynamo effect. Since the typical state of the magnetic 
T field is a collection of intense flux tubes, this effect is of limited 
interest. However, it may be important early in the evolution of the 
galactic magnetic field, and it will play a large role in numerical 
l simulations. Finally, we note that the formation of flux tubes is an 
essential ingredient in any successful dynamo model for stars or 
accreiion disks. 

^We will consider an idealized situation in which there exists a 
turbulent cascade with a scale L and a turbulent velocity, on the 
scale jof V T . We will assume that the magnetic field has an rms 
Alfv£n speed where V A ~ V T . We will also assume that the typical 

scale of curvature for the field lines is L. These assumptions are less 
restrictive ^an they may appear. If the turbulent cascade actually 
extends to larger length scales and higher velocities, then the mag- 
netic field is dynamically insignificant on these larger scales and we 
can still confine. our attention to scales of size L or smaller. If the 
magnetic field is In a shearing flow, surrounded by turbulence of its 
own creation, then the near equality of V T and V A is guaranteed, as 
well as the curvature of the magnetic field lines on the scale L. 

The field lines wilftend to stretch at a rate ~V-j/L. If the plasma 
is highly conducting then the same amount of matter will be en- 
trained on a progressively longer and longer flux tube. In a stationary 
state this stretching will be balanced by the pinching off of closed 
loops. These loops will have aradius ~L and a longitudinal compres- 
sive force ~ pV^/L. This tension will be opposed, usually by turbu- 
lent stretching with a force of ~V|/L. Some large fraction of the time 
the loops will collapse. Regardless whether the internal pressure of 
the loop is dominated by the magnetic field or gas pressure the 
magnetic tension will decrease more slowly than the turbulent 
stretching force and the loop will collapse to a plasmoid ball, whose 
energy is slowly lost to microscopic dissipation. This process will 
tend to remove matter from the flux tubes at a rate of -Y-j/L, which 
is rapid and will produce largely evacuated flux tubes under almost 
any circumstances. If we start from a uniform or nearly uniform 
field, this process will end when the same amount of flux is divided 
into some number of intense flux tubes with a magnetic pressure 
equal to the ambient pressure and a local P of order unity or less. The 
final rms Alfvdn velocity will be the geometric mean between its 
initial value and the local sound speed. This increase will occur at 
a rate comparable to Vy/L, in agreement with the results of numeri- 
cal experiments [2,3]. 

What will be the typical radius of the individual flux tubes? A 
single flux tube with an internal Alfvdn speed of % ~ C s , and 
exposed to an ambient turbulent velocity of V T , will remain coupled 
to the fluid provided that r, < (Vj/V At ) 2 L. On the other hand, these 
tubes will impede the flow, and thereby reduce the ambient fluid 
velocity below V T , if the total number N is large enough that Nr t /L 
is greater than l. The requirement that the magnetic energy be 
divided intoN flux tubes is just the requirement that Nr J VJ,~ V jJL 2 , 
which implies that the flux tubes will not impede the flow if r, is 
comparable to, or greater than, L(V A /V Al ) 2 . We conclude that the 


